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ABSTRACT 

The origin and progenitors of short-hard gamma-ray bursts remain a puzzle 
and a highly debated topic. The most popular model to interpret short-hard 
bursts is the merger of two compact stellar objects (NS-NS or NS-BH). A closer 
analysis by Zhang et al. reveals that a small sample of GRBs that have the most 
robust evidence of the compact star merger origin (Type I) is not a fair represen- 
tation of the BATSE short-hard GRB sample. This raises the question whether 
all short-hard GRBs are Type I, and whether there is a strong contamination of 
GRBs due to massive star core collapses (Type II) in the short-hard GRB sample. 
We utilize a Monte Carlo approach to determine whether a merger progenitor 
model can self-consistently account for all the observations of short-hard GRBs, 
including a sample with redshift measurements in the Swift era (2;-known sample) 
and the CGRO/BATSE sample. We apply various merger time delay distribu- 
tions invoked in compact star merger models to derive the redshift distributions 
of Type I GRBs, and then constrain the unknown luminosity function of Type I 
GRBs using the observed luminosity- redshift (L — z) distributions of the z- known 
sample and the peak flux distribution (log — log P distribution) of the BATSE 
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sample. We find that it is difficult to reconcile a merger scenario progenitor 
with all the data. In order to satisfy both observational constraints, the allowed 
merger delay time scale is essentially zero, which suggests that most short GRBs 
should follow the star forming history of the universe. We explore the possibility 
that the observed short-hard GRB sample is a mix of Type I (which requires a 
merger time delay to define the ^-distribution) and Type II (whose ^-distribution 
tracks the star forming rate) GRBs. We find that the required Type II fraction 
is very high, up to 90%. We conclude that Type I GRBs only account for a small 
fraction of short-hard GRBs and that the majority of short-hard GRBs track the 
star forming history of the universe and are likely of Type II origin. 

Subject headings: 



1. Introduction 



Short-hard gamma ray bursts (GRBs) have been an enigma since the identification of a 
bimodal distribution in the CRGO/BATSE data by Kouveliotou et al. (1993). This study 
showed that GRBs are distributed into two populations with shorter bursts having a harder 
spectrum and long bursts a softer spectrum, leading to the short-hard/long-soft classifica- 
tion. This purely observational division has a fair amount of scatter and does not necessarily 
indicate the nature of the intrinsic progenitor of a burst. Progress in understanding the 
progenitors of both long and short GRBs was made following the discoveries of their re- 
spective afterglows and host galaxies. While long bursts have been rnore securely shown 



to be associated with the collapses of m assive stars (iHjorth et al. 20031 : IStanek et al. 2003 



Campana et al. 20061 : iPian et al. 20061 ). the identification of a progenitor type for short- 
hard bursts has not been as successful. The most popular model is a merger event be- 
tween two com pact stellar objects, be it two neutron stars (NS-NS) or a NS and a black 



hole (NS-BH) flLattimer fc Schramm 19761 : Eichler et al. 19891 : iNaravan et al. 19921 ). This 



is supported by observational evi dence of a lack of a supernova component in a handful o f 



short-hard GRBs to deep limits flHiorth et al. 20051 : ICovino et al. 20061 : iKann et al. 20081 1 



as well as the very important discovery of a handful of short bursts identified in non- 



ine ealaxies. such as GRB 050709 (IGehre 


Is et al. 2005: 


Bloom et al. 2006: 


Fox et al. 2005: 


Villasenor et al. 2005: 


Hiorth et al. 2005: 


Barthelmv et al. 2005: 


Bereer et al. 2005). Such 



observational evidence, however, is not ubiquitous for all short GRBs. In fact, most short- 
hard GRBs d iscovered later are found in star forming galaxies or to have missing hosts 
jBerger 2009h . 
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Prompted by the discovery of GRB 060614, a nearby long GRB without an associated su- 
pernova but with many properties consistent with a merger-type p rogenitor (iGehrels et al. 2006 
Gal- Yam et ah 20061 : livnbo et aL 2o"o6l : IPella Valle et al. 2006h . Zhang et al. (2007) sug- 
gested that the long vs. short classification of GRBs does not necessarily match the physical 
origins of massive star core collapses and mergers of compact stellar objects, respectively. 
They suggested naming GRBs with massive star and compact stellar object merger origins 
as "Type II" and "Type I" , respectively, so as to be differentiated from the traditional "long" 
and "short" terminology. This viewpoint was elaborated in Zhang et al. (2009). In that pa- 
per, Zhang et al. (2009) showed not only that some long GRBs can be of a Type I origin, but 
also that a good fraction (not only the Gaussian tail of the long population) of short GRBs 
could be of a Type II origin. They argued that the two recently discovered high-z GRBs 
with intrinsic short durations, GRBs 080913 and 090423, are most likely Type II bursts, and 
further suggested that some high-z, high-L short-hard GRBs can be of the Type II origin as 
well. The goal of this work is to investigate through statistical methods whether the data 
are consistent with the hypothesis that "all short/hard GRBs are Type I, i.e. related to 
compact star mergers". 

Many of the specifics of Type I bursts are loosely constrained. Among others, two im- 
portant properties are the form of the luminosity function and the distribution of the merger 
time delay scale r, which is defined as the time elapsed between the formation of the two 
main sequence stars (i.e. the epoch of star formation) and the coalescence between the two 
evolved compact stellar objects (N S-NS or NS-BH). Several studies haye endeavored to add 



constraints to thes e distributions (lAndo 20041 : iGuetta &: Piran 20051 : iGal-Yam et al. 2005 



Nakar et al. 2006bl : iGuetta fc Piran 2006h . For the merger delay time scale, usually a long 
delay is invoked, in the form of either a roughly constant delay (anywhere from 1-6 Gyr) 
or a distribution that is proportional to a power 7 of the delay time scale r. Nakar, Gal- 
Yam & Fox (2006b) constrained the delay distribution of merger events to r > 4 Gyr or 
a distribution oc t~^'^ or shallower, while Guetta and Piran (2006) concluded that this 
distribution can be modeled by a logarithmic delay or one with a constant delay, gen- 
erally on the or der of a few Gyr. Later, some short GRBs with much higher redshifts 
were identified ( iBerger 20071 : iGraham et al. 20091 ). which po sed a challenge to the mod- 
els invoking a long m e rger delay time sca l e. Other studies (IBelczvnski &: Kalogera 2001 



Belczvnski et al. 2OO2I : llvanova et al. 20031 : IPewi fc Pols 20031 : IBelczvnski et al. 2006h sug- 
gest that merger timescales should not only be concentrated to long "classical" timescales 
(IBhattacharya fc van den Heuvel 19911 ) but also include a prompt merger channel. These 
arguments stem from the details of the binary evolution process. Belczynski and Kalogera 
(2001) as well as Ivanova et al. (2003) and Dewi and Pols (2003) proposed a scenario where 
ultra-compact orbits can be achieved by an extra mass transfer event in the evolution of the 
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binary, further reducing the orbital size of the final system that produces the bursts. This 
can lead to explosions after on the order of 10s of Myr, as seen in their population synthesis 
analyses. This scenario is consistent with the fact that short GRBs are seen in both early 
type galaxies and s tar forming galaxies (IBelczynski et al. 20061 : IZheng fc Ramirez-Ruiz 20071 : 
Zhang et al. 20091 ). however, it is unclear whether the model can reproduce all the available 
data of short-hard GRBs. 

The luminosity function of Type I GRBs is even more sparsely constrained. Nakar et 
al. (2006b) assumed a simple powerlaw luminosity function and found that an index -2 
can fit the available data by the end of 2005. Guetta & Piran (2006) introduced a broken 
powerlaw, with the indices in low luminosity ~ —0.5 and in high luminosity ranging from 
-1 to -2. Both works also included the caveat that the observational sample is very small, 
and that the data allows for some flexibility when combining rates, luminosities, and delay 
distributions. Modifications can also be added by considerations of multiple populations of 
short-hard bursts, such as a dual-peak luminosity fun ction to account for local SGR giant 
flare events (ITanvir et al. 20051: IChapman et al. 2009i), or the contri butions of Type I GRBs 
from globular clusters (iGrindlay et al. 20061 : ISalvaterra et al. 20081 ). but these contamina- 
tions are either not significant (Nakar et al. 2006a) or without robust evidence for their 
existence. 



Since the early at tempts of constraining the Type I progenitor models (INakar et al. 2006b 
Guetta fc Piran 20061 ) shortly after the discovery of the short GRB afterglows, the sample 
of short-hard GRBs with redshift measurements has significantly expanded. The observa- 
tions after 2005 suggest that the nearby, early type galaxies are not common short GRB 
hosts, and that a significant fract ion of short GRBs are likely from the high redshift universe 
( iBerger et al. 20071 : iBerger 20091 ). Furthermore, a closer analysis by Zhang et al. (2009) 
reveals that the Type I Gold Sample, a small group of GRBs that carry direct evidence in 
favor of a compact stellar object origin, are not representative of the BATSE short-hard 
GRBs. In particular, 4 out of 5 GRBs in the sample have extended emission, and all five 
have moderate hardness ratio. Even when not accounting for the extended emission, the 
"short spike" of the GRB 050724 has a duration longer than 2 seconds. This GRB is the 
"smoking gun" burst to establish the compact star merger origin of short GRBs, since it s 
afterglow was within a nearby early type galaxy (IBarthelmy et al. 20051 : IBerger et al. 20051 ) . 
However, after 4 years of observations, this burst is still the only one with such a robust 
associatior0. Zhang et al. (2009) suspect that some (maybe even most) short-hard GRBs 



^GRB 050509B |Gehrels et al. 2003 : iBloom et al. 20061 ) is believed to be associated with a cD elliptical 
galaxy in a nearby cluster. However, the argument was based on a chance coincidence argument with the 
XRT error box, since no optical afterglow was detected for this burst. 
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are not related to compact star mergers (Type I) , but are related to deaths of massive stars 
(Type II). Recently, Levesque et al. (2009) reported a study on z = 2.6 short-hard GRB 
090426 and concluded that the data are more consistent with the burst being of a Type II 
origin. Another related study is Nysewander et al. (2009) (see also Kann et al. 2007, Kann 
et al. 2008) which analyzed the optical afterglow properties of both long and short GRBs 
and found that the two populations have very similar optical-to-X-ray flux ratios, suggesting 
that the average circumburst density of the two populations is similar. This lends support 
to the conjecture that a good fraction of short-hard GRBs may not be of the Type I origin. 

The current sample of short-hard GRBs with redshift measurements (the 2;-known sam- 
ple) is large enough to serve the purpose of constraining redshift distribution and luminosity 
function of the population. This sample, together with the BATSE short-hard GRB sample, 
can be used to perform a self-consistent check of the hypothesis that "all short-hard GRBs 
have compact binary merger progenitors (Type I)" . This is the goal of this paper. Assum- 
ing that all the known short/hard GRBs detected in the Swift era are Type I GRBs, we 
perform a scries of Monte Carlo simulations to constrain the luminosity function of Type 
I GRBs by adopting several merger time delay distribution models, and check whether the 
same model can reproduce the peak flux distribution (log — logP) of the BATSE sample. 
In §2 we detail the model assumptions and some information on the simulations. In §3 we 
present the results of the constraints on various models with different combinations of merger 
time delay distribution and luminosity function. Our results are summarized in §4 with a 
short discussion. 



2. Models and Theoretical Framework 

An advantage of utilizing numerical methods to approach a problem is that it can be 
easily broken down into its constituent parts for easy processing. The number of any type 
of GRBs that occur within a comoving volume element, dV/dz, per unit observed time at 
redshift z z + dz and luminosity L L + dL is given by 

dN _RGRBiz)dVizl^^^^^ 



dtdzdL 1 + z dz 

with the factor of {1 + z) accounting for the cosmological time dilation, Rgrb{z) being the 
GRB volume event rate (in unit of Gpc~^ yr~^) as a function of z, ^{L) the luminosity 
function, and dV{z)/dz the comoving volume element given by 



dV{z) c AttD 



dz Ho{l + zy[nM{l + z)^ + QAY/^' 



(2) 
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for a flat A cold dark matter (ACDM) universe. Throughout the work, Q^n and Qa are set 
to 0.3 and 0.7, respectively. Both the expressions for the GRB volume event rate {Rgrb{z)) 
and the luminosity function {<^{L)) of Type I bursts are unknown, and are the major focus 
of this study. 

First, there is no theoretical prediction on the form of the luminosity function of Type I 
GRBs, which should in principle depend on the distributions of the masses of each member 
in the merging binary, the impact parameter, the collimation angle, and the viewing angle. 
Taking other astrophysical objects (e.g. long or Type 11 GRBs) as examples, we assume two 
general forms for the luminosity function for Type 1 GRBs, i.e. a powerlaw (PL) 



where a, cti, and a2 are the power law indices, L;, the break luminosity, and $o a normal- 
ization constant. 

Second, the ^-distribution of Type I GRBs can be modeled theoretically. For Type 
II GRBs, it is usually assumed that the volume rate of GRBs follows the star forming 
history, since the delay between the formation and death of a massive star is on the order 
of a few million years, much shorter than the variations in the cosmic star forming history 
or cosmological timescales. For Type 1 bursts, however, there is a delay, r, between the 
creation of the binary system, which follows the star forming history, and the eventual 
decay of the binary's orbit via gravitational radiation until coalescence. So the redshift 
distribution of Type I GRBs can be determined by the star forming history distribution 
convolved by a distribution of the merger delay time scale r. This latter distribution is 
not fully estabhshed theoretically, and we test the following three models that have been 
discussed in the literature: 

• Constant delay with dispersion: A 5-function hke delay of 1-5 Gry after the creation 
of the burst with the addition of a small normal dispersion of a — 0.3, providing for 
small variations in binary systems. 




(3) 



or a smoothed broken powerlaw (BPL) 





• Logarithmic delay: Delay where the distribution of delays P(log(r))dlog(T) ~ const, 
which implies P(t) ~ l/r. This form has been proposed by Piran (1992), Guetta and 
Piran (2006), and Nakar, Gal- Yam and Fox (2006). 
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Delay distribution from population synthesis: Belczynski et al. (2008) have modeled 
the NS-NS and NS-BH merger delay time scales using their population synthesis code. 
We have used their data and fit it with a 5*^ order polynomial and used this em- 
pirical model in our simulations. This model includes the fast merger channels and 
allows for many short (< 100 M yr) merg ers (see Fig{T]). We do not consider evolu- 
tionary effects due to metallicity (ILi 20081) or increase in the event rate with redshift 
( iKistler et al. 20071 : iKistler et al. 20081 ). since the former effect may not be relevant for 
compact star merger progenitors, and since Type I events are often more local events 
so that the redshift evolution effect, if any, is not pronounced. 



As will be evident later, none of the above models can interpret all the short GRB data. We 
are then forced to consider the possibility that some or even all short GRBs are not Type I 
events but are instead related to massive stars (Type II). We therefore consider the following 
two redshift distribution models as well. 



• No delay: All short-hard events have no delay from the star forming history, implying 
a non-merger (Type II) origin. 

• Mixed Type I/II distribution: The observed short GRBs are a combination of a Type 
I population (with a delay distribution defined by the population synthesis model, 
Belczynski et al. 2008) and a Type II population that follows the star forming history. 
The fraction of bursts in each population is a free parameter and can be diagnosed 
from the data. 



Once the value of the delay is assigned (in units of Myr) it needs to be added to the 
previously simulated redshift information to determine the redshift of the actual merger 
event, and consequently the redshift of the Type I GRB. Using the simulated redshift for the 
creation of the binary system, -Zcreation, which is assumed to follow the star forming history 
(i.e. SF2 model of Porciani and Madau, 2001), the cosmological look-back time is calculated 
by 

r ^creation C ^creation ^ "J^ 

With the values of this integral discretized over the simulated range (z=0-10) in units of 
Myr, we then subtract the merger time delay and re-convert to a redshift value with the 
same table. Those bursts with a negative lookback time (i.e. those that have not occurred 
yet) are discarded. The new redshift serves as the redshift of the merger and the GRB, -Zgrb- 
Figure [2] shows how the redshift distribution is affected by different models of the merger 
timescale distribution. 



- 8 - 



With this formahsm as a backdrop, we have enough information to create a set of bursts, 
each one defined with a unique and random (L, z) pair, which can then be passed through a 
series of filters that mimic a detector and then compared to the observed distribution. These 
simulations are similar to those conducted in Virgili et al. (2009), and here we summarize 
the most significant points. 

Fundamentally, Monte Carlo simulations rely on random numbers, and although it is 
possible to create true random numbers with a device that uses a stochastic process (e.g. 
thermal noise), pseudo-random numbers are much more convenient in terms of ease of use and 
possibility for exact repetition of simulations. In this code we utilize the SIMD-oriented Fast 



Merssene Twister, created by Mutsuo Saito and Makoto Matsumoto (ISaito fc Matsumoto 20081 ) 
of Hiroshima University. It is specifically designed for use with scientific Monte Carlo simu- 
lations, producing long strings of random numbers with a period of anywhere from 2®°^ — 1 

to 2216091 _ 1. 

In order to compare the simulated output with observations, it is necessary to be in 
the same band as the detector with which one is observing (i.e. the fc-co rrection). In orde r 



to achieve this, we assume every burst has a Band Function spectrum (IBand et al. 19931 ). 
with spectral indices a = —1.0 and (3 = —2.3 below and above a characteristic energy Eq. 
Although some BATSE and Swift bursts have been adequately fit with an exponential cut- 
off power law (Ghirlanda et al. 2004; Sakamoto et al. 2008), this is very likely because the 
flux above Ep is too low to constrain the high energy photon index of the Band function. 
Recent Fermi observations suggest that for bursts bright enough, the Band function can fit 
the spectrum for both long (e.g. 080916C, Abdo et al. 2009a) and short bursts (e.g. GRB 
090510, Abdo et al. 2009b). In the case of GRB 080916C, the Band spectrum extends 
several orders of magnitude in energy and supports our choice of the intrinsic spectrum 
for the simulated bursts. The characteristic energy correlates with the peak of the vF{v) 
spectrum by the relation Ep = Eq{2 + a) , which we assign from the relation proposed by 
Liang et el (2004) @ 

Ep/200keV = C(L/10^Vg s"^)^/^ (6) 

where C is randomly distributed in [0.1,1]. This energy can then be used for the fc-correction 
from the simulated bolometric luminosity in the rest- frame 1 — 10^ keV band (based on a 
certain luminosity function) into an arbitrary detector bandpass spanning the energy range 



^Although this correlati on was derived for long GRBs, it has b een found that short GRBs share the 
similar Ep — L correlation ( Ghirlanda et al. 2009t Zhang et al. 2009 ). The Ep — E^^iso correlations for the 
two categories are very different, mainly due to the shorter durations of short GRBs with respect to long 
GRBs. 



- 9 - 



(ei, 62). The fc-correction parameter is defined by 

•ioV(n 
1/(1+^) 



jll EN{E)dE 



(7) 



The last step is to incorporate the detector threshold condition. For CGRO/BATSE, we take 
a count rate trigger with a moderate sensitivity ~ 10"'' erg cm~^ s~^ and use it to derive the 
log — log P distributions of various models. Considering the imaging trigger capability of 
Swift/BAT, one needs to convert flux into fiuence in order to apply the detector sensitivity. 
In order to assign a random duration based on a T90 distributio n for short bursts we us e 



the BATSE population derived value of Tg^"""* = 0.33 ± 0.21 s (IKouvehotou et al. 19931 ). 
where 0.33 is the center of a distribution with gaussian scatter. The details of modeling the 
detector threshold are summarized in more detail in Virgili et al. (2009). 



3. Results 

The set of short bursts with known redshift is, to date, relatively small. We utilize the 
samples presented in Zhang et al. (2009) as our guide to attain a sample of short/hard GRBs 
as candidate Type I bursts, i.e. a "Type I gold" sample that includes 5 GRBs and an "other 
short/hard" sample that includes 20 GRBs. The "Type I Gold" sample bursts are those that 
either (1) have a late type or elliptical host galaxy or (2) have a localization with a large 
offset from the center of its host galaxy with relatively low specific star formation rate as well 
as stringent upper limits on the an associated supernova. The "other short/hard" sample 
include other short/hard GRBs detected in the Swift era. Since we focus on those short 
GRBs with secure redshift measurements, we remove seven bursts (GRBs 000607, 051210, 
060313, 060502B, 061201, 070809, and 080503) from the list, which either do not have derived 
luminosities or the claimed redshift is very uncertain (e.g. 051210). This leaves a total of 20 
bursts, whose properties are summarized in Table 1 and displayed graphically in Figl31 We 
assume that the redshift values of this sample are all correct, but c aution about the sma ll 



chance of mis-identification due to afterglow/host chance coincidence (jCobb fc Bailyn 20081 ). 
Among the highest redshift GRBs in this sample, GRB 070714B has z = 0.923 (Graham 
et al. 2009), and GRB 090426 has z = 2.6 (Levesque et al. 2009). GRB 060121 has two 
uncertain redshifts z = 1.7,4.6 (de Ugarte Postigo et al. 2006), and we take the smaller 
value z = 1.7. Some studies (e.g. Berger et al. 2007; Berger 2009) suggest that there are 
more short bursts at these high redshifts. If they indeed exist, as elaborated below they only 
strengthens the argument presented in this paper, namely, not all short GRBs are of the 
Type I origin, and that most short GRBs track the star forming history of the universe. 
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Utilizing the models presented in the previous section, we ran various sets of simula- 
tions combining different luminosity functions and merger time delay distributions. Varying 
the luminosity function parameters for a particular merger timescale distribution, one can 
compare the model predictions of GRBs in the L — z space with the L — z distributions of the 
data (see Fig|3]). Each test utilizes the ID KS probability in z {Pks,z), ID KS probability 
in L {Pks,l), and the total KS probability, defined as PKS,t = Pks,z x Pks,l- For a single 
power law LP (Eq.[3]), varying a leads to a distribution of PKs,t and the maximum PKs,t 
defines the most likely a value. Since the observed short GRBs seem to have an upper cutoff 
in the L distribution, it is more reasonable to adopt a smoothed broken power LF (Eq.[4]). 
However, since there are three free parameters (ai, ^2, and Lh), it is difficult to constrain 
all three parameters in detail. We therefore fix 0:2 = 2.5, and constrain ai and Lt together 
using the L — z criterion. This results in a series of PKS,t contours in the (ai, L5) plane, as 
shown in FigsS] and from which we can infer the best fit parameters for that particular 
model. Using these parameters, we can then construct a simulated 2D L — z graph with 
both the observed and simulated data plotted. 

Our next constraint for the simulated bursts is consistency with the BATSE logA^ — 
l ogP distribution. The sample is of 309 bursts selected from the BATSE 4B catalog 



(iPaciesas et al. 19991 ) with T90 < 2s (on a 64 ms timescale). When binned, the distribution 
gives a power law of slope of -1.12, extending up to about 1-50 ph cm~^ s~^, disregarding 
the turnover at low photon flux, which is an artifact of the detector. For any ^-distribution 
model, we adopt the most probable luminosity function derived by the L — z constraint, and 
simulate the log — log P distribution in the BATSE 50-300 keV band. The simulated pho- 
ton flux output was screened by a simple cut at the BATSE threshold of 0.25 ph cm~^ s~^. 
The model results are compared with the observed data. This allows for a self-consistent 
check of every model through the joint L — z and log A^ — logP analysis. Two examples of 
such self-consistency checks are presented in Figs. [7] and El 

Next we break down the results by merger time delay model and comment on the 
constraints imposed by the observations. 



3.1. Constant merger time delay 

We start with the constant merger delay models (with small scatter, Fig.2b). The five 
models we considered are r = 1,2,3,4,5 Gyr. These models are already disfavored by the 
host galaxy data of short GRBs: most short GRBs have star-forming host galaxies similar 
to those of long GRBs (Zhang et al. 2009 and references therein). However, since these 
models are simple, they allow us to study in detail how the general consistency with the 
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data depends on r. As shown in Fig.4 (a.-c, for r = 1, 3, 5 Gyr), all values of constant r can 
show some consistency in L—z^ implying the pre-break luminosity function slope of anywhere 
from ai ~ 0.3 — 1.2. The general trend is that ai increases (steepens) with the shortening 
of r. Except for the 1 Gyr model, all other models give a very shallow ai. Combining with 
the corresponding 2;-distribution, this always leads to a very shallow log — log P that is 
inconsistent with the data. The reason is that the merger models give a clustering of bursts 
at low- 2; (Fig.2), so that the shape of the luminosity function carries significant weight in 
defining the shape of log — log P. This is in contrast to the case of Type II GRBs, which 
are spread in a wide range of z so that both the redshift distribution and luminosity function 
play an important role to define log A^ — log P %. The r = 1 Gyr model is not ruled out by 
the L — z criterion, and can reproduce the log A^ — log P well. However, as stated above, this 
model is in contradiction with the host galaxy data. 



3.2. Logarithmic and Population synthesis 

More realistic Type I models are those that introduce a distribution of the merger delay 
time scales. We discuss the L ogarithmic distribution model and a more detailed model 



involving population synthesis (IBelczynski et al. 20061 ). These models have the advantage 



to allow for a range of merger times that can lead to observations of bursts in diverse host 
galaxy types. Both models affect the redshift distribution in similar ways and have similar 
results in all tests, and are therefore discussed together. The constraints from the observed 
L and z distributions all imply a very shallow luminosity function, generally have a slope of 
-0.2 or even larger, with moderate consistency (40-50%) for a broken power law luminosity 
function (see Fig|7j). The reason for a shallow luminosity function is to avoid overproducing 
nearby low-L Type I GRBs. These shallow luminosity functions severely overproduce at 
the high photon fluxes, giving a much shallower log A^ — log P curve than observed. These 
models are therefore not favored by the short GRB data. 



3.3. No delay 

The above analysis leads us to conclude that the hypothesis that "all short GRBs 
detected by BATSE and Swift are of the compact star merger (Type I) origin" is not justified, 
and that one needs to seriously explore alternative models. The first possibility is that all 



•^In an extreme, for a pure Euclidean geometry, the log N — log P always has a slope of -3/2 regardless of 
the shape of luminosity funciton. 
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short GRBs are actually related to massive stars and track the star forming history of 
universe. We apply the same observational constraints to this "no delay" model. The 
implied slope of the luminosity function is, as expected from the general trend of steepening 
with decreasing delay time, steeper than the lowest constant r model, weighing in at ai = 
1.42. The consistency with the observed L — z distributions is low, about 20%, but not 
sufficient to completely rule out an association. The implied logA^ — logP has a slope 
steeper than observation. A breakdown of the various tests is shown in FigJHl which suggests 
that attributing all short GRBs to Type II events is not justified. This is expected since 
several Type I Gold Sample GRBs have been identified (Zhang et al. 2009), suggesting that 
at least some short GRBs should be of the Type I (compact star merger) origin. 



3.4. Mixed population model 

We are forced to consider the possibility that the observed short GRBs include both 
compact star merger (Type I) events and events that are associated with massive stars (Type 
II). Although this possibility looks ad hoc at first sight, it is already implied by the data. 
Zhang et al. (2009) discussed various criteria to judge the progenitor systems of GRBs and 
concluded that Tgo is not necessarily an informative parameter to define the physical category 
of a GRB. After discussing a series of multiple observational criteria, they applied the criteria 
that more directly carry the progenitor information to define a Type I Gold Sample. They 
found that the Gold Sample bursts are relative long (and most have extended emission) and 
not particularly hard, and that they are not a fair representation of the short hard GRB 
sample. On the other hand, none of the high luminosity short/hard GRBs have been found 
in elliptical or early type galaxies. Instead they all fall in star forming galaxies. Zhang et al. 
(2009) therefore speculated that some or even most high-luminosity short GRBs are of the 
Type II origin. Our above analysis strengthens that conclusion and demands a more serious 
investigation of the mixed population in short GRBs. 

We test a merger timescale distribution that is the combination of a population of Type 
I bursts that have a range of distribution of the merger delay time scales as predicted by 
the population synthesis model of Belczynski et al. (2008), and a second component that 
tracks the star forming history of the universe (Type II). For the second component, we fix 
its luminosity function as the one inferred from the known Type II population, as discussed 
in Virgili et al. (2009) and Liang et al. (2007). The logarithmic distribution is very similar 
to the population-synthesis derived distribution, so we take the results from the latter as 
similarly applicable to the former. 

We start with the most extreme case with 100% Type II GRBs. This is similar to the 
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"no delay" model except that the luminosity function is fixed as that determined from the 
Type II GRB analysis (Virgili et al. 2009; Liang et al. 2007) rather than derived from the 
L — z constraint. This model shows a very good consistency with the logA^ — logP data 
(FiglHb), fit by a power-law with slope of -1.11 @ This case, however, is securely ruled out 
by the L — z analysis, with the probability contours for the luminosity function parameters 
showing no consistency between the simulated results and the observed sample, even with 
the inclusion of the recent high redshift bursts. Increasing the fraction of mixing with Type 
I GRBs results in a shallower log A^ — log P distribution (Figinb), but also a general increase 
in the consistency with the L — z constraints (FigjS]). In order to understand the origin of a 
shallower log — log P with an increasing Type I fraction, we plot in Fig|9] the distribution 
of photon fiux with merger time delay distribution. One can see the large grouping of bursts 
from the Type II population with zero time delay (solid dots) as well as a distribution of 
bursts with varying values of r contributing to all levels of photon flux. Most importantly, 
the bursts contributing to the very high photon flux end (> 100 ph cm~^ s~^), which is 
causing the log — log P to be shallow, are generally the bursts with merger time delays 
greater than 6 Gyr. The progenitors of these bursts are created very early in the universe but 
the GRBs do not occur until reaching a very local redshift, causing the high photon fluxes. 
In order to suppress these high flux bursts, one needs to keep a large fraction of Type II 
GRBs but only introduce a small contamination of Type I GRBs in the short GRB sample. 
Focusing on the L — z constraint alone, the best PKS,t is achieved at a moderate mix fraction 
between the Type II and Type I GRBs, reaching a 50% KS signiflcance for as much as 75% 
contribution of the Type II sample (FigJSt). However, to compromise the requirement of 
log N — log P flts, the most possible break down of the mix would be over 90% Type II and 
less than 10% Type I. For example, a contamination of 5% Type I bursts gives a 3% (4%) 
excess of bursts above 56 ph cm~^ s^^ (10 ph cm^^ s^^); the former being the highest photon 
flux bin in the BATSE log A^ — log P (see Figs J6] and [TOj) . A 10% contamination, which gives 
a smaller but non-negligible PKs,t, has an excess of 5% (7%) at the same benchmark photon 
fluxes. A breakdown of various constraints for the two mixed models (90-10 and 95-5) are 
shown in FigJTOl 



4. Conclusions and Discussion 

With the extensive afterglow follow up observations of short GRBs in the Swift era, a 
sample of z-known GRBs is signiflcantly expanded. This sample, together with the BATSE 



''There is a small excess at the high photon flux end. In view of the log scale involved, this only gives an 
excess of bursts above 100 ph cm~^ s^^ that is below 2. 
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short GRB sample, can be used to constrain the luminosity function and redsliift distribution 
of short GRBs by means of reproducing the observed L — z distribution and log — log P 
distribution. In this paper, we have performed a detailed analysis on a list of models using 
Monte Carlo simulations. Our results can be summarized as follows: 



The hypothesis that "all short hard GRBs are of a compact star merger origin (Type 
I)" is disfavored by the data. In particular, the popula r merger time delay model 



derived from population synthesis (IBelczynski et al. 20061 ) along with the logarithmic 



distribution model both demand a very shallow luminosity function in order to satisfy 
the L — z constraint. This is because a steeper luminosity function would over-produce 
low- 2;, low-L Type I GRBs that are not observed. Such a shallow luminosity function, 
combined with the redshift distribution derived from the merger delay time distribu- 
tion, leads to a very shallow predicted log N — log P, which is absolutely inconsistent 
with the BATSE logiV — logP data. Considering more unrealistic constant delay mod- 
els, those with large delays (2-5 Gyr) suffer the same problem. Only the 1 Gyr model 
can satisfy both constraints. However, this model predicts that most short GRBs re- 
side in early type galaxies with low star forming rate. This is ruled out by the host 
galaxy data of short GRBs, which shows a dominant star forming galaxy population. 

• A model that invokes no merger delay distribution can satisfy both the L — z and 
log N — log P better. This hints that most short GRBs may be related to massive stars 
(Type II). This model however predicts too steep a log — log P to be fully consistent 
with the data. 

• The data seem to demand a mix of Type I and Type II GRBs in the observed short 
GRB population. This is in consistent with the argument presented in Zhang et al. 
(2009) who argued that the Type I Gold Sample is not a fair representation of the 
BATSE short/hard GRBs. In order to get a steep enough logA^ — logP and in the 
mean time retain a non-negligible consistency with the L — z data, the required Type 
II "contamination" in the short GRB sample is very high, up to 90%, or even higher. 
The Type I GRBs (those associated with compact star mergers) only contribute to a 
small fraction (as low as 10% or less) of the short GRB population. 

We note that although the best mixed Type I/II model can marginally satisfy both the 
L — z and log — log P constraint, the significance of the criteria are not high enough for us 
to draw firmer conclusions. The hybrid model has too many parameters in the luminosity 
functions of the two components so that one cannot fully constrain them from the available 
data. In the analysis, we have adopted the luminosity function of known high-luminosity 
long GRBs (Virgili et al. 2009; Liang et al. 2007) as the one for the Type II component. 
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It is possible that the "short" Type II population has a different luminosity function from 
the "long" Type II population. This may allow a better consistency with the data. Further 
studies are needed to explore such a possibility. 

Our results imply that the rate of compact star mergers that give rise to short hard 
GRBs is smaller by near one order of magnitude than previously est imated based on th e 



assumption that all short GRBs are due to compact star mergers (e.g. iNakar et al. 2006bl ). 
The rate of gravitational wave bursts that are associated with compact object merger events 
is also lower by near one order of magnitude. 

Our results also demand that within bursts of the massive star origin (Type II), there is 
another factor that lead to the apparent bimodal distribution in the T90— hardness space. It 
may be related to the property of the progenitor star, or be related the process of launching 
a relativistic jet. The duration of a GRB is related to the duration of an relativistic jet tha t 



dissipates, which can be shorter than the total time scale of accretion (jZhang et al. 20091 ). 
Detailed models for short-duration Type II GRBs are called for. 

As shown, important and robust conclusions can be drawn about the nature of short/hard 
bursts from the current observations. Increasing the sample of short /hard bursts, especially 
those with redshift measurements and clear host galaxy associations, is of the greatest im- 
portance toward understanding the diverse underlying progenitors of these bursts and how 
we come to observe them. 
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Table 1: The 2;-known short/hard GRB sample 



GRB 


z 




name 


redshift 
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0.65 ±0.09 


060505 


0.0889 


~ 0.0009^ 


060801 


1.131 


47.6+f^ 
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Note, 
surements: 



uminosities deri ved by author unle ss otherwise specified 



GRB050509B: Gehrelseta: 



(jBarthelmv et al. 
GRB000607: 



2005 



GRB050709 
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References for redshift mea- 
(Fox et al. 2005h: GRB050724 



GRB050813: 
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GRB051210 



GRB 060121: jLevan et al. 20061 ) 



B loom et al. 2007l):GRB O6O5O5 

( Stratta et al. 2007 . ) 



GRB061210: 
GRB070714 B 
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(iBerger 20071): GR B061217: (jBerger 2007); GRB070429B: (|Cenko et al. 200: 

GRB 070724A: (Icucchiara et al. 20071): GRB071227 
( Levesaue et al. 2009h). GRB0905 10: |Rau e t al. 2009|): " We chose 



dCenko et al. 20081): 



z=1.7 for this analysis; Derived from 



on: (jPalmer et al. 2006h . E^^iso- (,Hullinger et al. 20061 ): 
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T (Myr) T (Myr) 

Fig. 1. — Fit to the Belczynski population synthesis results for the merger delay time scale 
(r) distribution of NS-NS and NS-BH binary systems. The left panel shows the fit to the 
data, while the right shows our Monte Carlo simulation output in gray superimposed with 
the population synthesis data in black. 



SFH 




Fig. 2. — Modified redshift distributions after taking into account the shifting from the 
merger timescale distribution, including the effects of cosmological time dilation and the 
comoving volume element, dV/dz. The left panel shows the unaltered star forming history 
(black), population synthesis (dark gray) and logarithmic models (light gray), while the right 
shows various constant delays. All histograms contain the same number of bursts. 
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Fig. 3. — Sample of the 2;-known short-hard GRBs detected in the Swift era. The "Type I 
gold" sample are shown as squares, while triangles are the "other short-hard" bursts. 
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Fig. 4. — A series of contours displaying the total KS probability, PKS,t for a sample of 
merger models and a model that has no delay from the star forming history. Panels (a)-(c) 
correspond to a 5 Gyr, 3 Gyr and 1 Gyr, (d) population synthesis model, (e) logarithmic 
model and (f) the no delay model. 
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Fig. 5. — A series of contours displaying the total KS probability, PKS,t, for various amounts 
of contamination with Type II in Type I GRBs. Left to right, top to bottom, (a)-(e): 20, 
50, 75, 85, 90, 95% Type II bursts, respectively. 



- 25 - 




1000 




P (ph cm V) [50-300 keV] 



Fig. 6. — Comparison of the log A^ — log P distributions for the various models. All, with the 
exception of a model with no delay from the star forming history and one composed entirely 
of Type II GRBs, are too shallow to accommodate the BATSE observations. The last pane 
(c) shows the population synthesis (dash) and logarithmic models (dot). The mixed models 
indicate how much contamination from Type II bursts are incorporated. The mix of 75%, 
for example, indicates 75% of bursts are chosen from the Type II luminosity function and no 
delay from the star formation history, while 25% are chosen from the population synthesis 
delay model. This allows us to constrain how contaminated the observed sample might be 
from Type II bursts. A 95% Type II mix shows a slope similar to the observed BATSE 
sample albeit with a small excess of bursts at the high photon flux end. 
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Fig. 7. — PKS,t contours, implied 2T) L — z distribution and log N — log P for the population 
synthesis (a & c) and logarithmic merger (b & d) time delay distribution models, together 
with a broken power law luminosity function. For the last pane (e) , the population synthesis 
distribution is the dashed line, the logarithmic the dotted. 
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Fig. 8. — A control simulation to show that a distribution of short bursts that strictly 
follows the star forming history (implicitly infer a Type II origin). The first three panes 
are the Pks,z, Pks,l, PKS,t contours. The last two are the 2D L — z diagram and BATSE 
band log N — log P implied from the luminosity function parameters at the maximum of the 
PKS,t contour. In the last panel, the observed BATSE distribution is the sohd fine while the 
simulated results are the dashed line. 
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Fig. 9. — Photon flux vs merger time delay for a simulation run with a broken power law lumi- 
nosity function and population synthesis merger timescale distribution contaminated by 80% 
high-luminosity Type II bursts, which follow the star forming history. The BATSE cutoff, 
which determines which bursts contribute to the log N — log P, is roughly 0.25 ph cm~^ s"^. 
Most of the bursts at high photon fluxes, which cause the log N-log P to be more shallow 
than observed, also correspond to the tail of the merger timescale distribution, having delays 
greater than ~ 5 Gyr. 
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Fig. 10. — A breakdown of the L~z distribution, PKS,t contour and log A^— log P distribution 
for 90-10 (a-b) and 95-5 (c-d) mix of Type II bursts and Type I bursts following a population 
synthesis derived merger timescale distribution. In the final panel (e), the BATSE sample is 
the solid, 90-10 mix the dashed, and 95-5 mix the dotted lines, respectively. 



